Complementarity of direct and indirect Dark Matter detection experiments 
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We investigate the prospects for reconstructing the mass, spin-independent and spin-dependent 
cross-sections of Dark Matter particles with a combination of a future direct detection experiments 
such as XEN0N1T, and the IceCube neutrino telescope in the 86-string configuration including 
the DeepCore array. We quantify the degree of complementarity between the two experiments by 
adopting realistic values for their exposure, energy threshold and resolution. Starting from bench- 
mark models arising from a supersymmetric model with 25 free parameters, we show that despite 
the stringent constraints set by the run with 79 strings, IceCube can help break the degeneracies 
in the Dark Matter cross-section parameter space, even in the unfortunate case where it fails to 
discover high energy neutrinos from the Sun. We also discuss how the uncertainties associated with 
the nuclear structure function and with astrophysical quantities such as the Dark Matter density 
and velocity distribution affect the reconstruction of the Dark Matter particle parameters from the 
combined datasets. 



I. INTRODUCTION 

The identification of the nature of Dark Matter (DM) 
particles is one of the central unresolved questions in 
modern particle physics and cosmology pQ. Weakly in- 
teracting massive particles (WIMPs) are among the best 
motivated candidates for Cold Dark Matter, as they are 
thermally produced and naturally lead to the correct relic 
abundance in the early Universe, and they naturally arise 
from well-motivated extensions of the Standard Model 
(SM) [M5]. 

DM can be searched for with a variety of detection 
strategies. In particular, WIMPs can be detected di- 
rectly, i.e. through their scattering off nuclei in un- 
derground detectors [5J Ej. Possible hints of detection 
of a light WIMP have emerged from data obtained by 
the CDMS II (SI, CRESST-II [§], CoGeNT QH] and 
DAMA/LIBRA [IT], but the interpretation of these 
events as due to spin-independent (SI) scattering of a 
10 GeV DM particle has been challenged [HI [J3] by sev- 
eral other experiments such as XENON100 [H], CRESST 
commissioning run [15] and CDMS [T6] [TT] . The next 
generation of direct detection experiments, characterized 
by very low background, is awaited to confirm or defini- 
tively rule out these ambiguous hints of detection. 

Even if new particles are convincingly discovered with 
direct detection experiments, reconstructing the proper- 
ties of DM particles, such as the mass and the WIMP- 
nucleon cross-sections, will be a complex task. First, the 
astrophysical parameters describing the DM distribution 
in the Galactic halo and in the solar neighborhood are 
affected by large uncertainties. It is known for instance 
that the DM velocity distribution may substantially de- 
viate from the idealized Maxwell-Boltzman (MB) distri- 
bution which is usually assumed [T5H20J . affecting the 
WIMP parameters reconstruction [3T] [52]. A realistic 
reconstruction of the DM parameters requires marginal- 
izing over astrophysical uncertainties [23-25], or inte- 
grating them out |26| . 

Second, the spectrum of nuclear recoils is insensitive 
to the WIMP mass when the latter is much larger than 



that of the nuclei of the target material, making a mass 
determination impossible for WIMPs heavier than ap- 
proximately 100 GeV (e.g. [27] [28]). Third, as recently 
pointed out in Ref. [29]) n ot only it is impossible to dis- 
entangle spin-independent (SI) and spin-dependent (SD) 
couplings with a single direct detection experiment, but 
the large uncertainties associated with the nuclear struc- 
ture function might lead to an error of about one order 
of magnitude on the reconstructed SD coupling [2l?l [3D] . 

Uncertainties and degeneracies in the parameter space 
can fortunately be reduced with a careful combination 
with other DM searches, e.g. combining the information 
arising from direct detection experiments with different 
targets [35] [5TH53] , or combining direct detection with 
accelerator searches [34]. Here we explore the degree of 
complementarity between DM direct searches and an in- 
direct detection strategy that consists of searching for 
high energy neutrinos produced by the capture and an- 
nihilation of WIMPs in the Sun's core with the IceCube 
neutrino telescope at the South Pole [35] (see also Refs. 
[361 137) for a similar analysis in the framework of super- 
symmetric models). 

In this paper, we quantify the degree of complemen- 
tarity of direct and indirect DM searches with a com- 
bination of a future direct detection experiments such 
as XENON1T, and the IceCube neutrino telescope in 
the 86-string configuration including the DeepCore array. 
We start by assessing the reconstruction capabilities of 
XENON1T [35] for 3 benchmark DM candidates within 
reach during the next 5 years: XENON1T is expected 
to reach 10 _47 cm 2 in sensitivity for the SI interaction, 
and 10 _42 cm 2 for the SD interaction (the same sensitiv- 
ity to SD as is expected for IceCube in the same time 
period 39J). We will demonstrate that even if one of the 
detectors does not see a signal the reconstruction of the 
physical parameters is still improved by utilizing both 
experiments. We will then assess the impact of uncer- 
tainties from astrophysics and from the nuclear structure 
functions in the reconstruction of the WIMP parameters. 
Our final results will be obtained by marginalizing over 
all these nuisance parameters, assessing the degradation 



of the reconstruction of physical properties, such as the 
DM mass and the WIMP-nuclei cross-sections with re- 
spect to the scenario with fixed nuisance parameters. 

The paper is organized as follows. We first introduce 
the theoretical expected signal for DM scattering off nu- 



clei in underground detectors, Section II A and for the 
neutrino flux from the Sun, Sec. Ill Bl Scctionllllldescribes 



the phenomenological approach we use in studying the 
WIMP signal as well as the statistical framework. Sec- 
tion llV] describes the set up of the future XENON1T ex- 



periment and its sensitivity for detecting WIMP signals. 
On the same lines, in Sec. [V] we describe the prospect 
for detection with IceCube and its sensitivity in recon- 



structing WIMP parameters. Section VI illustrates the 
effectiveness of combining different search strategies for 
reconstructing the fiducial benchmarks. We subsequently 
discuss the uncertainties that affect the signal reconstruc- 



tion in Sec. VII we first describe the impact of astrophys- 
ical uncertainties over the Galactic parameters, then the 
effect of undetermined nuclear structure functions, and 
lastly comment on the velocity distribution parametriza- 
Finally, we summarize our findings in Sec. 



tion. 



VIII 



II. PREDICTED SIGNALS FROM WIMPS 

In this section we review the theoretical predictions for 
the direct detection rate and for the neutrino flux arising 
from annihilation of DM particles in the Sun. The scope 
of this brief summary is to introduce the key model pa- 
rameters of our phenomenological analysis and underline 
the (different) dependence of the expect rates in direct 
and indirect detection experiments. 



A. Theoretical rate for direct detection 

Direct detection experiments aim to detect nuclear re- 
coils arising from the scattering of WIMPs off target nu- 
clei. The differential spectrum of a DM particle recoiling 
off a nucleus, in units of events per time per detector 
mass per energy, has the form 
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where E is the energy transferred during the collision, 
Po = pdm(-R©) is the WIMP density in the solar neigh- 
bourhood, mdm the WIMP mass, da/dE is the differ- 
ential cross-section for the scattering, and /(??(£)) is the 
normalized WIMP velocity distribution in the Earth's 
rest frame. The integration in the differential rate is per- 
formed over all incident particles capable of depositing 
a recoil energy of E. For elastic scalar interactions, this 
implies a lower integration limit of 



IMjsrE 
2/j, 
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where Mjj is the mass of the target nucleus, and [i = 
m-DM M_\f / (mnM + Ma/-) is the WIMP-nucleus reduced 
mass. As for the velocity distribution /(«(£)), we use 
the MB parametrization [101 HI] and neglect the time 
dependent modulation due to the Earth motion around 
the Sun. We defer to Sec. IVH Al and lVII Cl the discussion 
about the role of astrophysical uncertainties. 

The differential cross-section da/dE encodes the par- 
ticle and nuclear physics information and is in general 
separated into the SI and SD contributions as 
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a. Spin-independent interaction 
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where (x n — mjyM'rrin / {m~DM +Tn„) is the WIMP-nucleon 
reduced mass, af^ is the SI zero-momentum WIMP- 
nucleon cross-section, Z (A) is the atomic (mass) number 
of the target nucleus used, and f p (/„) is the WIMP ef- 
fective coherent coupling to the proton (neutron) . We as- 
sume the WIMP couples equally to the neutron and the 
proton (/„ = f p ), so that the differential cross-section 
da/dE is sensitive only to A 2 . The nuclear form factor 
J-si(E) characterises the loss of coherence for nonzero 
momentum transfer, and is well parametrised by the 
Helm form factor 021 BU for all nuclei HU 



Jsi(B) = 3e 

where s = 1 fm, r = \/R 2 



, sin(gr) — qrcos(qr) 
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"5s2, R = 1.2 A 1 / 3 fm, and 



q = y/2M^E. 

b. Spin- dependent interaction 
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where a^° is the SD zero-momentum WIMP-nucleon 
cross-section, a p (a n ) are axial WIMP-proton (neu- 
tron) couplings, J is the total spin of the nucleus and 
(S p ) ((S n )) is the proton (neutron) spin averaged over 
the nucleus. The nuclear form factor for SD is usually 
defined as 



Hd(E) = 



S(E) 
S(0) ' 



and 



S(q) = alS 0Q (q) + a aiS i(q) + a 2 Sn{q) . 



(7) 



(8) 



with ao — a n + a p (ai = a p — a n ) being the isoscalar 
(isovector) coupling. Furthermore we assume equal cou- 
pling to neutron and proton (a n — a p ), hence only the 
structure factor Sqo will be relevant for our analysis. 



This assumption is motivated by the fact that theoreti- 
cal model of WIMPs typically predict an similar cross- 
section to proton and neutron [531 E3 • 

For the xenon-based detector we will consider, there 
are two isotopes that have a non-zero total spin because 
of the unpaired neutrons: 129 Xe, with J = 1/2 and abun- 
dance 26.44%, and 131 Xe, with J = 3/2 and abundance 
21.18%. 

The structure functions Soo(q) and Su(q) are related 
to the transverse electric and longitudinal projections of 
the axial current. These functions can be computed in 
a shell-model for the atomic nucleus, and the spin of 
the nucleus is computed by means of the wave func- 
tions of the unpaired nucleons. Assuming a particular 
interaction between nucleons, these are placed in energy 
levels according to the exclusion principle and includ- 
ing as many excited levels as possible, making this kind 
of computation difficult. Finally the projected currents 
are computed by evaluating the matrix elements of the 
many-nucleon model. Nuclear shell models are more re- 
liable for heavy nuclei, but even in this case there can 
be significant deviations at zero momentum transfer or 
at high momentum. To bracket the uncertainties in the 
case of the xenon nucleus, for both isotopes we consider 
two parametrizations for the structure functions: the first 
(Nijmegenll hereafter) was computed in 1997 by Ressell 
and Dean [461 , while the second is based on a very recent 
computation using chiral effective field theory formalism 
and accounting for two body interactions [471 (from now 
on CEFT). In the first part of our discussion we assume 
nuclear structure functions from CEFT formalism, while 
in Sec. |VII B| we discuss the effect of marginalizing over 
various realization of the nuclear structure functions. 

The total number of recoils expected, as a function 
of the DM parameters, in a detector of mass Mdct in a 
given observed energy range [E m - ln , £ max ] over an expo- 
sure time T is obtained by integrating Eq. (nl) over energy 
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where we have accounted for an energy independent ef- 
ficiency factor e and a finite energy resolution u{E) for 
the detector: 
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For this analysis we will focus on the future XENON1T 
experiment, the details of which are outlined in Sec. |IV| 



B. Muon signal in neutrino telescopes 

Indirect detection of DM using neutrino telescopes in- 
volves four processes: capture of WIMPs by the Sun, 
annihilation of these WIMPs, production of neutrinos 
following the annihilation event, and finally detection of 



these neutrinos. The formalism described in this sec- 
tion is encoded in the DarkSUSY 5.0.6 software package 
which we use in our analysis [48] . 

Capture of WIMPs in the Sun occurs when WIMPs 
elastically scatter off nuclei in the Sun and lose enough 
energy to reduce their velocity to below the solar escape 
velocity. Subsequent scattering events reduce the velocity 
of the captured WIMPs further, until they concentrate 
and thermalize in the core of the Sun. The capture rate 
is H9] 



where 
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is the capture rate per unit shell volume, R is the solar 
radius and index i runs across nuclear species present in 
the Sun. The variable w is the velocity of the WIMP at 
the shell, and w = \/u 2 + v 2 , where u is the velocity at an 
infinite distance away from the shell {i.e. where the in- 
fluence of the shell's gravitational potential is negligible) 
and v is the escape velocity at the shell. The integration 
limit M ma x is given by 



TODM - Mjv; 



(13) 



The term wfl^^w) quantifies the probability that a 
WIMP will scatter to a velocity less than the escape ve- 
locity, and is proportional to ain-i, where tii is the number 
density of nuclei i in the shell and <j{ is the total interac- 
tion cross-section between the WIMP and nuclei i. 

The WIMP-nuclei cross-section Oi can be expressed 
as 0EO] 
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The Sun is composed overwhelmingly of spin- ^ hydro- 
gen and spin-0 helium, with only small quantities of heav- 
ier elements. The increased atomic number of these heav- 
ier elements will compensate for their lower abundance 
via the A 2 term and give appreciable contributions to the 
SI part of oi. However no such enhancement occurs for 
the SD part, with total spin being related to the number 
of unpaired nucleons. Thus we can consider only a SD 
contribution from hydrogen nuclei, reducing Eq. ( 14 1 to 
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s " for i = 1, 
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(16) 



While SI interactions are taken into account in the cap- 
ture rate calculation, in practice the SD interaction is 
dominant. As the SD process occurs directly on proton, 
the theoretical rate is not affected by nuclear structure 
functions describing the coherence of the nucleus. 

As shown in Eq. (11) the capture rate is dependent 
upon the density of WIMPs. However as WIMP capture 
is a continuous process it is sensitive not to the local den- 
sity at the Sun's current position but instead samples the 
local density along the prior path of the Sun around the 
galaxy [51 . For this analysis we assume a constant local 
DM density P0 bs , and we will discuss the uncertainties 
related to this assumption in Sec. |VH A| 

Acting against the accumulation of WIMPs is the pro- 
cess of DM annihilation. We can describe the total pop- 
ulation of WIMPs in the Sun N(t) [52] by the equation 



dN(t) 
dt 



= c c - r„(t) 



(17) 



The 



where T a (t) = ^C a N 2 (t) is the annihilation rate, 
parameter C a is dependent on the distribution of WIMPs 
in the Sun and (a a v), the zero velocity WIMP annihila- 
tion cross-section [2]. 



Solving Eq. (17) gives us an expression for the annihi- 
lation rate: 



r Q (i) = ^tanh 2 ^ 



(18) 



where r = 1/ \JC c C a is the capture-annihilation equilib- 
rium time scale. In the case of t 3> r, Eq. ( 18 1 reduces 
to T a (t) — \C C , and equilibrium between capture and 
annihilation is reached. For this study we assume this 
steady state scenario; it allows us to take a more model 
independent approach by eliminating the dependence on 
(cr a v)- 

WIMPs which have accumulated in the Sun can an- 
nihilate with each other and produce Standard Model 
particles. The majority of the decay products of these 
particles are absorbed almost immediately and without 
consequence in the core of the Sun. However, certain 
classes of WIMPs can decay directly into neutrinos, and 
in other cases the SM decay products can themselves de- 
cay into neutrinos, which can escape from the Sun and 
potentially be detected on Earth. Detection of these neu- 
trinos by neutrino telescopes occurs via the observation 
of the Cerenkov radiation emitted by the particles pro- 
duced following the weak force interactions between the 
neutrinos and the matter in and around the telescope. 
Of particular interest are muons created by the charged 
current interactions of muon neutrinos, as their range is 
such that they create long, relatively detectable tracks 
compared to other leptons [S3] . Optical detection of the 
Cerenkov radiation in a transparent medium such as wa- 
ter or ice then allows the incoming neutrino's energy and 
origin to be reconstructed. For this analysis we will focus 
on the IceCube neutrino telescope, the details of which 
are outlined in Sec. W\ 



III. STATISTICAL FRAMEWORK AND 
FIDUCIAL WIMP MODELS 

To illustrate the capabilities of reconstruction and 
complementarity of future direct detection experiments 
and neutrino telescopes we consider three fiducial bench- 
mark points described in Tab. [I] These WIMP models 
represent a phenomenological approach in the descrip- 
tion of the theoretical parameters O = {fflDM^nVn} 
we are interested in, and capture the relevant aspects of 
the analysis: 

• Benchmark A, is characterized by a light mass 
of 60 GeV and its cross-section on nucleus is dom- 
inated by the SD component. 

• Benchmark B has an intermediate DM mass of 
100 GeV and sizable WIMP-nuclei cross-section for 
both SI and SD interactions. 

• Benchmark C is characterized by a heavy mass of 
500 GeV and significant SI but negligible SD cross- 
sections. 

We have checked that our conclusions are robust and hold 
for benchmarks with same behavior of the cross-sections 
but different masses as well. These fiducial WIMP mod- 
els are representative of well- motivated neutralino config- 
urations arising in scans of the MSSM25 54 parameter 
space, which is a phenomenological parameterization of 
the minimally supersymmctric Standard Model (MSSM) 
with 25 free parameters defined at the electroweak scale. 

For each fiducial benchmark we simulate mock data, 
using the experimental likelihoods £(0), from the true 
model, i.e. without Poisson scatter [55]. These theoreti- 
cal signals expected in the detectors, which are the num- 
ber of recoiling nuclei in XENON1T arising from both 
SI and SD part and the number of up-going muons N^ 
in IceCube, are given in Tab. [I] and will be described in 
details in Sees. |IV| and [V] respectively. 

For the sampling of the theoretical parameter space we 
adopt the Bayesian methodology. We employ the public 
code MultiNest v2 . 12 [56, 57 , which uses an ellipsoidal 
and multimodal nested-sampling algorithm to map the 
posterior over the full parameter space: 



V{@\d) oc £(e)7r(6) 



(19) 



where 7r(0) is the prior probability density function 
(pdf). The priors for the three theoretical parameters 
moM, c^ 1 and <7^ D are chosen to be flat on logarithmic 
scale so as not to favour any particular order of magni- 
tude, and are defined as follows: 

log 10 (TO DM /GeV) : 1 ->• 3 , 
log 10 (aS7cm 2 ):-60^-43, 
log 10 (aS D /cm 2 ):-55^-38. 

We set n livo = 25000, an efficiency factor of 10" 4 and 
a tolerance factor of 0.01 156 , which ensure that the 



TABLE I. Below are listed the benchmark points used in the analysis, together with their mass, cross-sections, counts predicted 
in XENON1T and expected muon signal in IceCube (and the WIMP annihilation channel), as labelled. The first benchmark 
is characterized by a light DM mass and large SD contribution, the second has intermediate mass and both sizable SI and 
SD cross-sections, while the third has a large mass and dominant SI contribution. As for the SD contribution to XENON1T 
detector, we quote the predicted counts arising from both nuclear structure functions considered in this analysis (see text for 
further details). 
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sampling is accurate enough to have a parameter esti- 
mation similar to Markov-Chains Monte Carlo sampling 
methods. The resulting chains are analysed with an 
adapted version of the package GetDist, supplemented 
with MATLAB scripts from the package SuperBayeS [SSJ 
155] . Two-dimensional posterior pdfs, V ma ,ig, marginal- 
ized over the nuisance parameters and the remaining n— 2 
theoretical parameters, are obtained from the chains by 
dividing the relevant parameter subspace into bins and 
counting the number of samples per bin. An x% credi- 
ble interval or region containing x% of the total volume 
of Pmarg is then constructed by demanding that "P mar at 
any point inside the region be larger than at any point 
outside. The inferred pdfs are sensitive to the choice of 
the mass prior range, as we have checked enlarging our 
prior from 10 GeV up to 100 TeV: the x% contours suffer 
from volume effects related to the behavior of the likeli- 
hood at very large mass, above 10 TeV. We however argue 
that these effects are not relevant for our analysis, since 
for these high masses the DM most likely does not meet 
the WIMP requirements anymore. Hence the mass prior 
range used here is driven by the standard mass range for 
WIMPs. 

The astrophysical parameters describing the DM halo 
and the DM in the solar neighborhood are regarded as 
nuisance parameters. These parameters are the local 
standard at rest velocity vq, the DM escape velocity in 
the halo v esc and p© . All these quantities have uncertain- 
ties that range from 20% to a factor of two. Indeed the 
observed values with one standard deviation, from [60 
IB"?] , are given by 



„obs 

yObS 

^esc 
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230 ±24.4 km s , 
544 ±39 kms" 1 , 
0.4 ±0.2 GeV cm" 3 



In the first part of our analysis we keep these nuisance pa- 
rameters fixed to their observed value, while in Sec. | VII A"] 
we present the final results marginalizing over them. The 
uncertainty over the nuclear structure function for SD 
interactions is also treated as a nuisance parameter in 
Sec. IViTBl 



IV. RECONSTRUCTION WITH XENON1T 

XENON1T will be the ton scale continuation of the 
XENON100 detector, and will be constructed under- 
ground at the Laboratori Nazionali del Gran Sasso in 
Italy. Start of data taking is planned for 2015, with first 
results released in 2017. 

We consider the energy window for DM searches to be 
from 10 to 45 keV, divided into 7 bins. In this range, the 
XENON1T likelihood is given by the product of inde- 
pendent Poisson likelihoods describing the probability of 
observing the predicted events for a given WIMP signal 
over the energy bins labelled by the index i 
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XENON1T 
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(20) 

N° hs is the observed number of events in each bin 
and 9bg,i is the total background in each bin, de- 
scribed below. The detector energy resolution is <r{E) — 
0.6 keV y/ E /keV . The effective exposure e e g = 2 ton x 
year comes from e x 77 cut x Anr = 5x 0.8 x 0.5, where e is 
the total exposure, rj cu t is the cut efficiency and Anr is 
the nuclear recoil acceptance. As XENON1T is expected 
to be almost background free [55], we consider a back- 
ground of 4 x 10~ 8 counts/kev/day/kg, which accounts 
for a total of #bg,i = 0.02 events per bin after integrating 
over the exposure time. The background is taken to be 
energy independent without loss of generality 

The expected number of recoiling nuclei for our fidu- 
cial points is shown in Tab. [I] Benchmark A has a 
negligible SI cross-section of about 10 _49 cm 2 , which is 
well below the expected sensitivity of XENON1T of 
~ 10~ 47 cm 2 : the theoretical predictions are compatible 
with the background in this case. All the signal comes 
from the SD contribution, which can range from 200 to 
400 events, depending on the nuclear structure function 
adopted for describing the xenon nuclei. In this case 
ofj D = 2.0 x 10 _40 cm 2 is close to the upper bound of 
XENON100 for SD interaction on neutron [30] and in 
the sensitivity range of XENON1T for SD which is ex- 
pected to be around 10~ 42 cm 2 [38 . Benchmark B is 
similar to A as far as it concerns the SD contribution, 
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FIG. 1. The reconstruction with XENON1T of the A (top), B (middle) and C (bottom) benchmark points with fixed astro- 
physics. No nuclear uncertainties are taken into account. The left, central and right plot shows the 2D marginal posterior pdf 
in the {jtidm, 0"n I }-P lane ! { m DM, <r^ D }-plane and {o - ^ 1 , erjp }-plane respectively. The contours denote the 68% and 95% credible 
regions and the diamond point shows the fiducial values of the benchmark point. 



with a number of events ranging from 100 to 350 de- 
pending on the structure function model. In contrast to 
Benchmark A however, B has a large contribution of 252 
events from the SI interaction, which has a cross-section 
of af 7 = 8.8 x 10 _46 cm 2 , just below the exclusion bound 
of XENON100 P3]. The last fiducial point, Benchmark 
C, has a negligible contribution from the SD interaction, 
accounting for zero events, while the SI interaction pre- 
dicts about 70 events. 

Figure [l] shows the reconstruction of our fiducial points 
with XENON1T only, with the galactic parameters fixed 
at their observed values and with the CEFT nuclear 
structure function in use. Hence no nuisance parameters 



are included in this part of our analysis. The left panels 
show the 2D marginal posterior pdf in the {mDM,^'}" 
subspace, the central panels show the 2D marginal pos- 
terior pdf in the {mDM ) < T n D }" su k s P ace j while the right- 
most panels show the 2D marginal posterior pdf in the 
{ct„, cr„ }-subspace, for Benchmarks A, B and C in the 
first, second and third rows respectively. 

If we allowed the DM particle to have only one inter- 
action, SI or SD, contributing to the observed events, we 
would expect the reconstruction of the WIMP-nucleon 
cross-section to be effective, namely to select only a 
closed region in the parameter space. However, as soon 
as both cross-sections are allowed to vary the marginal 



posterior pdf exhibits a tail towards very small values, 
since each parameter can compensate for the other to 
achieve the observed number of events. In other words, 
as the signal in XENON1T receives contributions from 
both SI and SD interactions a degeneracy can occur. This 
is apparent from all plots, where the 95% credible region 
does not exhibit a closed contour. The effect is more pro- 
nounced for the SD interaction, while in the case of SI 
scattering the 68% credible region is denoted by a closed 
contour around the fiducial point for Benchmark A, or 
always bounded from below in the B and C cases. This 
is due to the fact that XENON1T is more sensitive to SI 
interactions than SD interactions. In addition a low sig- 
nal for one interaction will cause the corresponding true 
value to be poorly reconstructed: for example, for A the 
SI interaction is very small and close to zero observed 
events, hence the signal in XENON1T at 68% confidence 
level (CL) is attributed to the largest cross-section that 
accounts for such events, and the true value is contained 
only in the 95% CL. This effect is mirrored for Bench- 
mark C, where again the true value for the SD interaction 
is only contained in the 95% credible region, which are 
as well upper bound for that interaction. 

In the reconstruction of the DM mass parameter, we 
encounter the usual limitation of direct detection: after 
the WIMP mass matches and then exceeds the target 
nuclei mass, the mass reconstruction becomes progres- 
sively worse, as once past the nuclei mass the signal rate 
drops as 1/todm- The degeneracy in the mass recon- 
struction can be seen as we move to the higher mass B 
and C points, as denoted in the left and central panel 
of Fig. [T] by the tail extending towards large masses. 
The only exception is Benchmark A, whose posterior pdf 
does not exhibit this mass degeneracy and appears to 
be a closed contour around the fiducial point at both 
68% and 95% CLs. A viable solution to ameliorate the 
DM mass reconstruction is to combine different target 
materials, typically with different atomic number, as de- 
scribed in [251 E3] . However heavy WIMPs still exhibit 
a tail in the posterior pdf extending towards high mass 
values. We will show in Sec. VI that it is possible to 



break this degeneracy by complementing the signal seen 
in XENON1T with the signal of IceCube, or even with 
the non-detection of a signal. 

It is worth commenting on the features appearing in 
the marginal posterior pdf contours at the 95% CL, in 
particular affecting benchmark C of Fig. [I] These fea- 
tures are not physical, but are an indication that the 
posterior pdf is flat, which makes it a challenge to sam- 
ple effectively. This also occurs as with IceCube alone or 
when the likelihood is poorly constrained by the data. 



V. RECONSTRUCTION WITH ICECUBE 

The IceCube neutrino telescope consists of 5160 dig- 
ital optical modules (DOMs) arranged vertically along 
86 strings embedded in a cubic kilometer of extremely 



transparent natural ice below the South Pole [55] • In- 
cluded in IceCube is the DeepCore sub-array, which con- 
sists of six strings arranged at the centre of IceCube with 
closer horizontal spacing and instrumented with DOMs 
of higher sensitivity and closer vertical spacing along the 
string. DeepCore has been designed to reduce the energy 
threshold of IceCube from 100 GeV down to 10 GeV [53]. 

The WIMP annihilations can produce a wide vari- 
ety of final states, whose subsequent decays can produce 
high energy neutrinos, with the branching ratios to each 
state dependent on the properties of the DM particle 
and the underlying theory that generates it. In our phe- 
nomenological approach we consider three annihilation 
final states with branching ratio equal to one, as labelled 
in Tab. [ij. For Benchmark A we consider the t + t~ final 
state, which gives a harder neutrino signal than the bb 
channel. Benchmark B annihilates into W + W~ , which 
is the dominant channel for a common fermionic WIMP 
with a mass above the W threshold. Finally we consider 
the VfJ/p channel for Benchmark C. The annihilation di- 
rectly into neutrinos is heavily suppressed for Majorana 
particles {e.g. neutralinos) |70) but can occur for vector 
WIMPs, such as the Kaluza-Klein photon [31] ITT]. 

For each benchmark point and its annihilation channel 
we first calculate the theoretical signal rate using the de- 
fault routines in DarkSUSY 5.0.6 EH] [72] (which employ 



the methods outlined in Sec. II B ) with a slight modifica- 
tion to enforce the steady state scenario. The MSSM25 
models we build upon for our fiducial WIMP models, see 
Tab. HJ have equilibration times of the order of 1 — 4 x 10 8 
years, e.g. an order of magnitude smaller than the age 
of the Sun, making our assumption of a steady state sce- 
nario a valid one. We use an exposure time of five 180 day 
austral winter observing seasons {e.g. 900 days), and an 
angular cut around the solar position of cu t = 3°. The 
predicted muon signals S^(0) for each benchmark point 
and annihilation channel are given in Tab. [I] From the 
background files included in the DarkSUSY 5.0.6 release 
we extracted the number of background events within 
the 3° angular cut. For a 180 day season the number 
of background events was 41, yielding 205 background 
events for our 900 day exposure time. With this back- 
ground estimation we have an estimate of the number of 
events needed for a detection within 5 seasons of data 
taking in IceCube, which is approximately 50 events, de- 
pending on the annihilation channel. This is compatible 
with the analysis released recently by the IceCube col- 
laboration [35 : by analyzing the data from one winter 
season with up-going muons plus one summer season with 
down-going muons, taken using 79 strings (including the 
six DeepCore strings), the strongest exclusion bounds are 
for a WIMP mass of 250 GeV annihilating into gauge 
bosons and with a spin dependent cross-section on pro- 
ton of about 10 _40 cm 2 producing ~ 15 events. The use 
of DeepCore allows the limit to be extended down to 
a mass of 20 GeV but the sensitivity rapidly diminishes, 
and at this mass only a SD cross-section of approximately 
10 _38 cm 2 can be excluded. For a WIMP with a mass of 



50 GeV the current exclusion bound is close to 10~ 40 cm 2 , 
comparable with the bound obtained with XENON100. 
For our fiducial points the expected number of events 
range from ~ 67 muons for Benchmark B to 7.8 muons 
produced by the v^v^ final state for C. Benchmark A 
annihilating to the t + t~ final state is intermediate with 
22 expected muons. Comparing these to Tab. I of [35] we 



can see that all our benchmarks have signal rates com- 
patible with the current exclusion limits. 

The likelihood function we use for the IceCube ex- 
periment is presented in |73j and included in DarkSUSY 
5.0.6. The likelihood of observing 7V to t events given a 
background rate of (9bg and a theoretical model 0, which 
has a true expected signal S tTne {&) = eS^(Q) events, is 



£lC86(e) = £(Wtot|0BG,S,i) = 



27TOV 



JVto 



9 B G + £^(e)) JVtot e 



-(flBG+eS^e)) 
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-exp 



1 /lne 

2 1^7 



de. 



(21) 



The variable e quantifies potential systematic errors be- 
tween the true expected signal St r ue(©) and the nomi- 
nally predicted signal 5^(0). The relative fractional er- 
ror on 5^(0) is then e — 1, and this is then marginalized 
over in a semi-Bayesian way, assuming a log-normal prob- 
ability distribution and with standard deviation a e . This 
likelihood takes into account only the total number of 
signal events. A likelihood function taking into account 
the energy of the signal muons is also presented in |73j 
and encoded into DarkSUSY 5.0.6, but this requires an 
accurate description of the energy response function of 
IceCube, which is not publicly available for the 86 string 
configuration. 

There are several experimental aspects of IceCube that 
are salient to our discussions. First, the effective area of 
IceCube, which quantifies the sensitivity of the detector, 
increases with energy [71] ■ Thus higher mass WIMPs 
will tend to produce stronger signals as they will, in gen- 
eral, produce higher energy neutrinos and thus higher 
energy muons. Second, as neutrino energy increases the 
mean angular error for the reconstructed muon track de- 
creases [74] , During calculation of signal in IceCube an 
angular cut (</> cu t) around the solar position is made. 
Thus lower mass WIMPs will produce lower signal, as 
the muons they produce can potentially have greater an- 
gular deviation from the solar position and so fewer will 
survive the angular cut. Conversely the higher energy 
muons will be more densely clustered around the solar 
position, and so more will pass the angular cut, increas- 
ing signal for the higher mass WIMPs. 

Given these limitations in the number likelihood of Ice- 
Cube, there exists a degeneracy between WIMP mass 
and WIMP-nuclei cross-section tr t ot = 0"„ + <r^ D . Con- 
sider two scenarios, one with a WIMP of high mass and 
low cross-section ertot , and another with a WIMP of low 
mass but high Ctot- The first scenario will have a low 
capture rate, as this is dependent on er t ot and inversely 
dependent on WIMP mass, Eq. (11). Assuming capture- 



annihilation equilibrium, this scenario will feature a lower 
muon flux, but the muons will be of a higher energy and 
so hence produce higher signal due to the features of Ice- 
Cube outlined above. The muons of the second scenario 
will be of lower energy, but they will have a higher flux 



as their capture rate is higher. This can compensate the 
lower muon energy and generate a high signal. As the 
energy spectrum, which could be used to determine the 
DM mass, is not defined in the number likelihood, these 
two scenarios cannot be distinguished from each other. 
As for direct detection signals, there is a degeneracy be- 
tween er„ and c„ D . The solar capture rate, given in 



Eq. (11), receives contributions from both SI and SD in- 



teractions, however the sensitivity of IceCube to the two 
type of interactions is different, as shown in [35]: the SI 
is poorly constrained and the actual exclusion limits are 
one or two orders of magnitude less constraining than the 
XENON100 ones, while for the SD cross-section the ex- 
clusion limits are competitive with XENON100 at high 
masses. 

Fig. [2] shows the reconstruction of the fiducial model B 
with IceCube only, and exemplify the discussion above. 
In the left panel, which shows the 2D marginal posterior 
pdf in the {raDM^^j-plane, the signal clearly exhibits 
an almost independent behavior with respect to SI inter- 
action. On the contrary, the SD cross-section as a func- 
tion of the DM mass is both lower and upper bounded 
at 68% CL (middle panel) . The mass degeneracy is a re- 
sult of the missing spectral information. The degeneracy 
between cr„ and <r„ is similar to the case of direct de- 
tection, as shown in the right panel and in in Fig. [Tl As 
far as it concerns Benchmarks A and C, the muon sig- 
nal produced by WIMP annihilations in the Sun cannot 
be disentangled from the background, leading to a flat 
and featureless marginalized posterior pdf and underlin- 
ing that there is not detection in IceCube. Once again, 
the features in the pdfs have no physical meaning and 
are merely symptoms of a flat likelihood. 



VI. COMBINED XENON1T AND ICECUBE 

ANALYSIS 

Figure [3] shows how the future detection/constraint 
from IceCube complements a detection in XENON1T, 
for fixed astrophysical parameters. The top row of Fig. [3] 
illustrates the combined reconstruction of Benchmark 
A. Even though this benchmark would produce only 12 
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FIG. 2. Reconstruction using IceCube (denoted IC86), with fixed astrophysics, for benchmark point B with annihilation channel 
W + W~ . The left, central and right plot shows the 2D marginal posterior pdf in the {moM, fj„}-plaiie, {mDM,o"n D }-plane and 
{<x„ , er„}-plane respectively. The contours denote the 68% and 95% credible regions and the diamond point shows the fiducial 
values of the benchmark point. 



events in IceCube, hence too small a signal to claim a 
detection, it still has the ability to improve the recon- 
struction. Comparing to the top row of Fig. [T] the 68% 
and 95% CLs have shrunk. The reconstruction of er^ 1 and 
fj^ D exhibits the same trend as the case of XENON1T 
only, except that now the 68% CL for SD cross-section 
in the right panel is upper bounded, because of the Ice- 
Cube upper bound. The mass determination remain un- 
changed, demonstrating the good reconstruction capabil- 
ity of XENON1T and the reduced sensitivity of IceCube. 

The fiducial point B (second row of Fig. [3]) fully ex- 
ploits the properties of complementarity between DM 
search strategies and is the principal illustration of the 
main point of our analysis, as it can be detected by both 
experiments. Comparing with Figs. [I] and [2] we can see 
that the combination of the two experiments allows for 
the SD cross-section (middle panel) to be well determined 
at the 68% CL. The mass degeneracy is also significantly 
reduced and it is well determined at the 68% CL. The 
contours for the SI as a function of jtidm have been con- 
tracted sensibly with respect to the separate experiments, 
but at 68% CL the SI scattering cross-section is only 
upper-bounded, as can be seen as well from the third 
panel: the SD is determined within one order of magni- 
tude at most and contains the true value, while the SI 
can not be larger than 10 _44 cm 2 . Therefore if a WIMP 
is detected by both experiments, meaning that it has siz- 
able cross-sections, the SD cross-section and the mass can 
be well reconstructed, while a precise reconstruction of SI 
interaction would require the use of a further experiment; 
in this case a direct detection experiment with different 
target material would be the optimal choice [33] . 

Benchmark C is illustrated in the last row and shows 
remarkably the constraining power of IceCube, even with 
only 8 events from high energy neutrinos. The most 
dramatic improvement is in the determination of cr^ 1 : 
where previously the 68% and 95% CLs extended from 
the upper bound all the way down to the bottom of our 
prior range, they are now bounded from below as well. 



As expected, only an upper limit can be set on a™ : 
the 95% CL is clearly bounded from above at values of 



T SD 



a," 10 cm . The effect of the SD upper bound is much 
more striking than for Benchmark A 7 signaling the fact 
the IceCube is more likely to be sensitive to high WIMP 
masses. Note that in this analysis we assumed, in absence 
of published information on the energy response function 
of the 86-string configuration, that nothing can be said 
on the energy of the neutrinos measured by IceCube. It 
is known that the energy resolution of IceCube is poor; 
for a muon of energy 10 4 — 10 8 GeV in the core of the 
detector there is an uncertainty of a factor of two [75] . 
Nonetheless, including the spectral information on ob- 
served events would certainly improve the mass determi- 
nation, one of the main weaknesses of the reconstruction 
procedure from XenonlT data only, and substantially re- 
duce the degeneracy in the DM particle parameter space. 



VII. ROLE OF UNCERTAINTIES IN THE 
RECONSTRUCTION OF WIMP PARAMETERS 

In this section we first assess the effect in the com- 
bined reconstruction of fiducial points when DM galactic 
parameters vary within their uncertainties (see Sec. Ill) 
and discuss the e ffectiv e DM density probed by the cap- 
ture rate (see Sec. II B I. We then account as well for our 
lack of knowledge on the true form of the nuclear struc- 
ture function (see Sec. II A). At the end of this section, 



we return to the additional uncertainties induced by the 
shape and parametrization of the velocity distribution. 



A. Astrophysical Uncertainties 

A reconstruction of fiducial points fixing the astrophys- 
ical parameters to their observed value could hold for 
very well measured parameters, however, as mentioned 



in Sec. Ill these are known only up to a various degrees 
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FIG. 3. Reconstruction using XENON1T combined with IceCube (denoted IC86), with fixed astrophysics, for benchmark 
points A, B and C with annihilation channels t + t~ (top), W + W~ (centre) and v^v^ (bottom). The left, central and right 
plot shows the 2D marginal posterior pdf in the jmDM,fff} -plane, {moM, a;p}-plane and {a^ 1 , a^ D } -plane respectively. The 
contours denote the 68% and 95% credible regions and the diamond point shows the fiducial values of the benchmark point. 



of precision. We therefore introduce these uncertain- 
ties into our analysis, treating the galactic parameters 
as nuisances V'astro and marginalizing over them. The as- 
trophysical likelihood governing these parameters follows 
a Gaussian distribution centered on the observed value. 
The total likelihood is then 

C = £(V>astro) X £XEN0N1T(S, Castro) X £t C 86(@, V'astro) • 

(22) 
The results of our analysis are illustrated in Fig. |4j which 
shows the 68% CL for marginalized astrophysics denoted 
by the medium thickness red line, for A (top), B (cen- 
ter) and C (bottom) . The primary effect of marginaliz- 
ing over the astrophysical parameters is to broaden the 



contour regions, which can be clearly seen by comparing 
the red contours to their fixed astrophysics counterparts 
given by the thick light green line. 

There is an additional uncertainty related to the lo- 
cal DM density and to the different way this parame- 
ter enters into the IceCube and XENON1T rates. The 
direct detection rate depends on p Q at present time at 
the Sun's position, while the solar capture rate depends 
on the long term history of the Sun, which completes 
an orbit around the Galactic Center in ~ 2 x 10 8 year. 
If we drop the assumption that the DM galactic halo is 
isotropic and smooth, throughout its journey the Sun will 
cross over-dense or under-dense regions with respect to 



11 



an averaged density ((p©)), to which the capture rate is 
sensitive. This difference has been evaluated for a triaxial 
DM halo arising from N-body simulations [76J . In case of 
simulation with baryons the difference is of the order of 
30%, while it can reach a factor of 3 for pure DM simula- 
tions. The MSSM25 models which serve as the starting 
points for our fiducial benchmarks have small annihila- 
tion cross-section, close or below the thermal one, hence 



the equilibration time scales r (see Sec. II B ) ranges from 
9 x 10 7 year for benchmark B, up to 2 — 3 x 10 9 year 
for points A and C . These values imply that the cap- 
ture process is fully sensitive to this averaged DM den- 
sity (case B as well, since it spanned half of the solar 
period). Even if we do not know the initial position of 
the Sun we could assume the extreme scenario such that 
the local DM density is different from the average density 
and corresponds to a over-dense or under-dense region. 
In this case we would get a 30% or even a factor of 3 
bias in the theoretical predictions. In other words the 
signal in one detector could be potentially boosted or 
suppressed by the same factor with respect to the signal 
in the other detector, as both rates scale linearly with 
the local DM density, affecting the complementarity in 
either a positive or negative way. 



B. Nuclear Uncertainties 



The structure functions presented in Sec. |II A| are im- 
portant quantities to be computed. Regarding the xenon 
nuclei we are interested in, the variation of Soo(g) leads 
to a different prediction of the events in the detector and 
hence, for example, to a huge uncertainty in the exclu- 
sion bound [3D]. Considering our fiducial benchmarks, 
the difference in the counts for the two structure func- 
tion is illustrated in Tab. [T] The Nijmegenll structure 
function predicts fewer events in the detector with re- 
spect to the CEFT structure function: the variation is 
sizable for A and B, where the events are affected by a 
factor of 3 of uncertainty while it is completely negligi- 
ble for case C. The reconstruction in XENON1T of a 
benchmark point is therefore affected by the assumption 
of the structure function: if the reconstruction of <r^ D 
parameter is done assuming a different form factor with 
respect that used for the mock data a systematic offset is 
generated. If the reconstruction is done assuming CEFT 
and the mock data are generated with Nijmegenll, this 
leads to an underestimation of the SD cross-section to 
account for the largest number of observed events; vice 
versa the SD cross-section will be overestimated. 

To account for the nuclear uncertainties in our 
Bayesian framework we follow the procedure of [29], 
which bracket these uncertainties with a formula with 
three free parameters 



^o(«) = tfi((l-&)e- a ' tt + &), 



(23) 



where j = 129 or 131, depending on the xenon isotope, 
Nj is a overall normalization factor, olj drives the slope in 



the low momentum regime, while f3j controls the height 
of a possible tail at large momentum. The variation of 
all parameters encompasses the whole region contained 
between the CEFT and Nijmegenll parametrizations but 
it is not meant to be a proper fit of the nuclear structure 
functions. We take a flat prior for all these additional nui- 
sance parameters ip nnc \, over their allowed range, which 
is: 

JV 129 = 0.045 -> 0.070 and JV 131 = 0.025 -> 0.052 , 

ai29 = 3.8 -» 4.0 and a 131 = 3.8. -» 4.0 , 

0129 = 0.013 -> 0.029 and /3 131 = 0.10 -> 0.12 . 

We therefore generate for each benchmark point new 
mock data based on the mean value of these parameters 
and, as before, we assess the impact on the reconstruction 
of the fiducial benchmarks by letting them vary within 
their prior range, in combination with the astrophysical 
uncertainties. The new likelihood is given by 

C = £(V>Tot) X £xENONlT(©,V ; Tot) X £lC86(@, V'astro) , 

where i/'Tot accounts for both astrophysical and nuclear 
uncertainties. Notice that the latter do not affect the 
reconstruction in IceCube: in solar capture the dominant 
contribution for SD arises from the hydrogen atom, which 
does not need a coherence factor for describing its single 
proton. Thus IceCube potentially has the capability of 
resolving the nuclear uncertainties when combined with 
XENON1T. 

In Fig. [4] the results for the 2D posterior pdf marginal- 
ized over all nuisance parameters are denoted by the dark 
blue solid line at 68% CL. This contour is not shown for 
Benchmark C as the nuclear uncertainties are not con- 
sidered because of the negligible a^° cross-section. The 
main result is that the nuclear uncertainties are smaller 
or at most comparable with the astrophysical ones. This 
can be explained by the fact that we are considering the 
Soo{q) nuclear structure function, which is known to be 
the better determined nuclear structure function, hence 
minimizing the role of nuclear uncertainties. Our findings 
are compatible with the analysis illustrated in [5D1 [33] . 



C. Effect of the shape of the velocity distribution 

The shape of the DM velocity distribution affects the 
direct detection and IceCube signal rates in different 
ways. The nuclear recoil threshold for DM detection in 
XENON1T is in the range of 10 keV, which implies that 
only particles with the largest velocities have sufficient 
energy to produce a detectable nuclear recoils above, 
as given by the v m ; n definition in Eq. (2l). For a 100 
GeV WIMP, as we have for Benchmark B, the typi- 
cal u m i n to produce a nuclear recoil is of the order of 
^min ~ 300 km s _1 . This means that it is crucial to accu- 
rately describe the high velocity tail of the DM velocity 
distribution in the galactic halo. On the other hand, as 



explained in Sec. II B the solar capture rate is sensitive 
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FIG. 4. Reconstruction using XENON1T combined with IceCube (denoted IC86) for benchmark points A with annihilation 
channel t + t~ (top), B with annihilation channel W + W~ (middle), and C with annihilation channel v^Vp, (bottom). The 
left, central and right plot shows the 2D marginal posterior pdf in the {m_DM, cr^j-plane, {ttidm, ajpj-plane and {a^, <7^ D }- 
plane respectively. The contours denote the 68% credible region as labelled: the green line (light thick solid) stands for fixed 
astrophysics and fixed SD structure function (SDSF), the red curve (darker solid) denotes marginalized astrophysics, while the 
dark blue outer line stands for all nuisance marginalized over. The diamond point shows the parameter values of the benchmark 
point. 



to the low velocity tail of the DM velocity distribution, 
with the maximum velocity of a WIMP that can be cap- 
tured being given by Eq. Q. For a 100 GeV WIMP the 
maximum velocity is of the order of u max ~ 200 km s _1 . 
Thus it is important to describe the low tail of the veloc- 
ity distribution as well. This is another physical quantity, 
besides the DM density in our galaxy, biasing the direct 
detection and indirect detection rate in different ways. 

Several N-body simulations of Milky Way-like galaxies 
suggest that the velocity distribution deviates from the 



standard MB halo. To deal with the discrepancy in the 
high velocity tail for direct detection experiments it has 
been proposed in Ref. [77] a velocity distribution modeled 
as a function of one parameter k and based on the these 
simulations. 

In our analysis, to assess if the shape of the velocity 
distribution is relevant for our analysis, we follow a dif- 
ferent approach, which allows us to consider the effect of 
both tails of the velocity distribution. We consider the 
Aquarius N-body simulation [2fJ with only DM, and the 
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FIG. 5. Bias due to parametrization of DM velocity distribution, benchmark B with varied astrophysics. Each contour denotes 
the 68% CL for the fiducial data point generated with an N-body simulation or with the MB distribution, as labelled, and 
reconstructed assuming a standard halo model. 



simulation by Ling et al. [TB], based on the RAMSES 
code, which includes baryons as well. We fit the mean of 
the N-body simulation Aquarius with a MB distribution, 
keeping the escape velocity fixed to 565 km s" 1 as given 
by the simulation but letting ^o vary: the best fit point 
underestimates both tail of the velocity distribution. For 
Benchmark B the difference in the number of events is of 
the order of 17% in XENON1T and of 15% in IceCube. 
By allowing a similar procedure with v csc = 520 km s _1 
in the case of Ling et al. the velocity distribution is much 
closer to the MB shape, and the discrepancy in the num- 
ber of predicted events is only of 4%. The mean of these 
numerical velocity distributions is used to generate new 
mock data for our benchmarks in the case of XENON1T 
and rescaled for IceCube. These new fiducial points are 
then reconstructed assuming as usual a MB velocity dis- 
tribution. We applied the procedure for Benchmark B, 
which is the point that takes full advantage of the com- 
bination of XENON1T and IceCube. 



Figure [5] shows that our analysis is robust against vari- 
ation of the shape of the velocity distribution, when 
both experiments are combined. From left to right 
the 2D marginalized posterior pdf for Benchmark B is 
shown, in the {mDMiC^j-plane, {mjM, cr^ D }-plane and 
{cr^jCT^j-plane. The middle thickness blue curve is the 
68% CL for a reconstruction of the DM physical parame- 
ters assuming that the true velocity distribution is given 
by the Aquarius simulation: the discrepancy with the 
MB curve (thick light orange) is very limited and affects 
mostly the SI cross-section with a most stringent upper 
bound, closer to the true value. The discrepancy with the 
Ling et al. distribution is very small, as the 68% contour 
level (thin black line) follows closely the MB contour. In 
other words the shape of the DM velocity distribution 
will likely have a limited impact in the reconstruction of 
WIMP physical parameters, in case of a positive detec- 
tion. 



VIII. CONCLUSIONS 

We have discussed the reconstruction capabilities for 
the physical WIMP parameters with a combination of 
data from XENON1T and IceCube. We have focused our 
attention on the degeneracy between the reconstruction 
of the SD and SI contributions, and we have shown the 
degeneracy between these two contributions arising from 
the analysis of the two experiments separately can be 
lifted from combining the two data sets. 

To illustrate and quantify the degree of complementar- 
ity, we have focused on three benchmark models: Bench- 
mark A, characterized by a light DM mass and large SD 
contribution; Benchmark B has intermediate mass and 
both sizable SI and SD cross-sections, and Benchmark 
C, with a large mass and dominant SI contribution. 

The combined analysis for Benchmark A allows us to il- 
lustrate the fact that even in cases where the parameters 
are such that no signal can be observed with IceCube, 
taking this fact into account in the global likelihood does 
improve the reconstruction of the DM parameters, as can 
be seen from the top row of Fig. [T] The reconstruction of 
fj^ 1 and <7^ D in this case exhibits the same trend as the 
case of XENON1T only, except that now the 68% CL 
for SD cross-section in the right panel is upper bounded, 
because of the IceCube upper bound. The mass determi- 
nation remains unchanged, demonstrating the good re- 
construction capability of XENON1T and the reduced 
sensitivity of IceCube. 

The fiducial point B (second row of Fig. [3]) fully ex- 
ploits the properties of complementarity between DM 
search strategies and is the principal illustration of the 
main point of our analysis, as it can be detected by both 
experiments. Comparing with Figs. [T] and [2] we can see 
that the combination of the two experiments allows for 
the SD cross-section (middle panel) to be well determined 
at the 68% CL. The mass degeneracy is also significantly 
reduced and it is well determined at the 68% CL. The 
contours for the SI cross-section as a function of ttidm 
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have been contracted sensibly with respect to the sepa- 
rate experiments, but at 68% CL the SI scattering cross- 
section is only upper-bounded, as can be seen as well 
from the third panel: the SD is determined within one 
order of magnitude at most and contains the true value, 
while the SI can not be larger than 10 _44 cm 2 . Therefore 
if a WIMP is detected by both experiments, meaning 
that it has sizable cross-sections, the SD cross-section 
and the mass can be well reconstructed, while a precise 
reconstruction of SI interaction would require the use of a 
further experiment; in this case a direct detection experi- 
ment with different target material would be the optimal 
choice [55] . 

Benchmark C, finally, shows remarkably the constrain- 
ing power of IceCube, even with only 8 expected events 
from high energy neutrinos. The most dramatic improve- 
ment is in the determination of <7„: where previously the 
68% and 95% CLs extended from the upper bound all the 
way down to the bottom of our prior range, they are now 
bounded from below as well. As expected, only an upper 
limit can be set on er^ D : the 95% CL is clearly bounded 
from above at values of <r^ D 10~ 42 cm 2 . The effect of the 
SD upper bound is much more striking than for Bench- 
mark A, signaling the fact the IceCube is more likely 
to be sensitive to high WIMP masses. In the absence 
of a publicly available energy response function for the 
86-string configuration of IceCube our analysis could not 
utilize the energy information of the neutrinos observed 
by IceCube, limiting its mass determination capability. 
Including this spectral information on would certainly 



improve the mass determination, one of the main limi- 
tations of the XENONlT-only reconstruction, and thus 
substantially reduce the degeneracy in the DM particle 
parameter space. 

The inclusion of uncertainties on astrophysical param- 
eters and nuclear structure functions only slightly affects 
the quality of the reconstruction of the 68% CLs, the 
degradation being worst in the case of less constrained 
scenarios, such as benchmark C, while B appears to 
be fairly robust against variation of astro and nuclear 
physics parameters. In other words the combination of 
XENON1T and IceCube in case of a positive signal could 
also be a way of constraining the DM velocity distribu- 
tion in the galactic halo, once all the subtle effect of the 
astrophysical parameters enters into the capture and di- 
rect detection rate are consistently taken into account. 

The strategy of combining direct and indirect searches 
for WIMPs allows for a better reconstruction of WIMP 
parameters, but the identification of the DM candidate 
will remain a difficult task. To understand the nature of 
DM it is necessary to complement the direct or indirect 
search strategies with LHC searches, as discussed e.g. in 
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